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Abstract 

Three decay mechanisms are derived systematically from the QCD 
Lagrangian using the field correlator method. Resulting operators 
contain no arbitrary parameters and depend only on characteristics of 
field correlators known from lattice and analytic calculations. When 
compared to existing phenomenological models, parameters are in 
good agreement with the corresponding fitted values. 


1 Introduction 

An enormous amount of experimental data on strong decays of mesons and 
baryons is partly used by theoreticians for comparison in the framework of 
the 3 P 0 model (T[] , and its flux-tube modifications 0. The analysis done in 
0 confirmed the general validity of the model, whereas in (4j results of other 
forms of decay operators have been also investigated in meson decays, and in 
|| in baryon decays. On the whole, the phenomenological picture seems to 
be satisfactory for the 3 Po model with some exclusions discussed in |3j and 
0. The recent extensive study of strong decays of strange quarkonia based 
on the 3 P 0 model was done in J7J. 

The key element which is missing in this situation is the systematic deriva¬ 
tion of all terms in the decay Hamiltonian from the basic principles, i.e. from 
the QCD Lagrangian. It is the purpose of the present paper to make some 
progress in this direction using the Field Correlator Method (FCM) § and 
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background perturbation theory || to treat nonperturbative (NP) QCD con¬ 
tributions together with perturbative ones. 

In doing so one should take into account the special role of pions in the 
hadron decays and therefore to perform accurately the chiral bosonization 
of the effective quark Lagrangian, obtained from the basic QCD Lagrangian. 
This will give the first term in the decay Hamiltonian, and the corresponding 
decay mechanism will be referred to as a Chiral Decay Mechanism (CDM). 
At the same time one should take into account the string degrees of freedom 
in the original meson and the possibility of string breaking due to the qq 
pair creation. The corresponding term in the decay Hamiltonian will be 
derived below without free parameters and this second mechanism will be 
called the String Breaking Mechanism (SBM). As will be seen, the dominant 
term of SBM has the structure, which can be compared quantitatively with 
the phenomenological fits of the 3 Po model in |3|. [|. 

Finally, the QCD perturbation theory in the perturbative background 
developed in |J . allows to derive two additional terms in the decay Hamilto¬ 
nian: one for the OZI-allowed decays, which has the Lorentz form of the 3 Si 
type but proceeds through the intermediate hybrid state, and another for 
the OZI-forbidden decays, which proceeds through the intermediate glueball 
state and only at very small distances reduces to the two-gluon or three-gluon 
qq pair creation. We shall call these mechanisms the Hybrid Mediated Decay 
(HMD) and the Glueball Mediated Decay (GMD) respectively. 


2 The Chiral Bosonization and the Chiral De¬ 
cay Mechanism 

One starts with the QCD Lagrangian in the Euclidean space-time and aver¬ 
ages over the gluonic fields writing the general form of the gauge-invariant 
correlator (known from lattice or analytic calculations, see refs, in [§) which 
contains the confining part D(x), namely, 

jy-(tr(F^(x)$(x, y)F Xa (y)$(y, x))) = {8^ x 8 va - 8^8 vX )D(x - y) + 0{D Q 

( 1 ) 

where 0(D i) contains a relatively small nonconfining part D\(x) and $(£, y ) = 
Pexpig jy A^dz^ is the parallel transporter. 
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Assuming also that all higher correlators can be neglected (as it is sup¬ 
ported by lattice data, see 0) one obtains the effectice action 0 
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-( 2 ) 

J EQL 


(e/) = / d4x j rf V^a Q (^) / ^(^)(^) 9 ^(l/) 9 ^a £ (l/)x 


X 


[l*hie J E{x, V ) + V)} 

where the kernel J(x,y ) is expressed through D(x), 

Je(m)(x, y) = J ^ dui dviK^ M)) {u, v)D(u -v), i = 1, 2, 3. 

and D(x) is connected to the string tension a in the usual way, a = \ J d 2 xD(x). 
Here the string kernel Je(m) contains color electric (magnetic) helds and the 
former is dominant at large distances, and therefore magnetic part will be 
omitted for simplicity. At this point one should specify the choice of in¬ 
tegration contours, the kernel K and the point Y in As it is known, 
the full result of the integration over gluon helds does not depend on the 
gauge and on the shape of the contours (when all correlators are taken into 
account); however to write the contribution of bilocal correlator ([[]) in the 
gauge-invariant form, one has to use one of the variants of the contour gauge, 
e.g. [|T^1 , where = 1, and to choose the contour corresponding to the 
minimal string which minimizes the contribution of higher correlators. In 
this section we consider the following geometry: the quarks at the points x, y 


( 2 ) 


(3) 


in (§) are at one end of the string (they are dynamically close |]TT[), while 
the point Y is at the position of the heavy antiquark, as it was introduced in 


11[. In the next section we shall consider a more general geometry, when the 
points Y in (^) are different and integration over dui and di\ runs over two 
different pieces of the broken string. For the results of the present section 
the exact definition of the point Y is inessential since the pion is emitted 
from the end of the string under consideration while another end of string is 
a spectator. 

As the next step the bosonization of the Lagrangian (|2|) can be done in 
the usual way, however with nonlocal mesonic helds with the result [O] 


A L = i / d A xd 4 y'ijj + (x)M(x,y)'ilj(y) 


(4) 


where the kernel M ( x , y) can be written as a nonlinear form for PS meson 
helds 

M(x,y)=Ms(x,y)exp(i'y 5 }(x,y)) + ..., 0=--^-, F n = 93MeV. (5) 
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Here ellipsis implies all other omitted terms, including those containing 
isovector scalars, vector and pseudovectors. 

One should note, that in our case when confinement is present the con¬ 
stant condensate of scalar-isoscalar field always enters multiplied by J(x,y) 
and thus produces the scalar confining potential of the string in M$(x,y ) 
Namely, for long enough string, i.e. 
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for lx — Y| 


gluonic correlation length in D(x), one has approximately |lT 

M s (x,y) w a|x|h (4) (x -y), 


T g , where T g is the 


( 6 ) 


This is different from the instanton, or the NJL model, where the scalar- 
isoscalar field acquires a nonzero condensate, which is constant in all space- 
time. 

Inserting (6j) into (^) one obtains a localized form of the quark-meson 
Lagrangian, describing interaction on one end of the string; to the lowest 
order in PS field, 


A\A 


A = [ ^(x)cr|x — Y| 75 ^— il>(x)dtd 3 x, 

J iv 


(7) 


One can visualize in 0 the simultaneous presence of quark fields ip, ip, 
together with the string cr|x — Y| and Nambu-Goldstone (NG) fields (j) A , 
A — 1,..., n 2 f — 1. 

Using Dirac equation for the quark field ip(x) one arrives as in [|T^] at the 
familiar Weinberg Lagrangian [|14[ 


A L c = g q A tr{ip^^ b uj^ip),uj g = ~{ud g u + - 


( 8 ) 


where u(x) = exp(|q 5 (p(x,x)). 
the local limit, 


In our derivation g q A is uniquely defined in 


9\ = 1 


(9) 


which agrees with large N c limit, discussed in ||T4 . 

Note that both Lagrangians 0.® are local limits of nonlocal expression 
0), and for not very long string with the length L ~ T g ~ 0.2 fm the 
nonlocality is essential. However string is apparently not present in (||), both 
quark operators there are solutions of Dirac equation with the string entering 
as a scalar potential. It is clear that the Lagrangian (|) describes the pion 
field emission both from the quark at one end of the string and from the 
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antiquark at another end of the string. In case of baryons one should sum 
up in (H) over all three quarks. 

The form (|p was used in |H| for the calculation of pionic transitions in 
the heavy-light mesons with g q A playing role of a fitting parameter, which 
turned out to be around 0,7. This difference from (|9j) can be considered as 
an indication of a possible role of nonlocality. 


3 The String Breaking Mechanism 


This mechanism was considered in some detail in [|16|] (see also refs, therein). 
In the present section we shall consider the pair creation vertex due to the 
nonperturbative QCD configurations. In the 3 P 0 model [l|] this vertex was 
modelled by an adjustable constant and in ]16[ by some function F. It is our 
purpose here to derive this vertex from the basic Aq effective action (|5J). 

To describe the creation of the qq pair in the presence of the string which 
connects quark Q and antiquark Q, one can as in [|TT[ take the large N c limit 
of the same Aq Lagrangian (H), which obtains by replacing 1-0^ (rc) 9 ^^?/) —> 
5fgN c Sp^y(^X, y). 


The resulting Lagrangian describes the creation of the qq pair at the 
points x, y respectively. Physically it is clear that x and y should lie on 
the (deformed) string connecting quark Q at the point X and antiquark 
Q at the point X. Since the final two strings connect X and x, and X, 
and y, we can also choose the contours of integration along the strings: e.g. 
A^x) = f* duiEi(ui), A^y) = J|. dviEi(vi), and the string is along the 
axis 1. (This change of contours from (||) can be traced to the effect of 
cancellation in the sum of contours from the quark Q to the point Y and 
with opposite sign -from the antiquark q to the point Y , which results in the 
contour integral between positions Q(X) and q{x)). As a result the kernel 
J E is 


<^e(x, x 4 ; y, yA) = / du x f di^D^i - rq) = -exp(- ^ 4 ) (10) 


ix 
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where for D(x) the Gaussian form was used, c.f. [|lTJ . As in (|7]) one can find 
the effective mass operator (due to color electric fields) M(x,y), using the 
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M(x,y ) = -i'uS{x,y)') A J{x,y) 


( 11 ) 
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and the estimate of the quark Green’s function S(x,y ) for a long string, 
L T g , done in ]Tl], [L3|, gives S(x,y ) iS^ 3 \x — y). 

As a result one obtains for the effective Lagrangian the same form as in 


A L < ' SBAI} —if d 4: xd 4 y^ + (x)AI < ' br \x,y)^ a (y) 


where 


M (br \x ) = -<5 ( 3 ) (x-y)e" 


(*4-S4r 
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Integrating over d(x 4 — y 4 ) one gets in the Minkowskian space-time 
al (sbm) = ^ f d ± x ^ x ^ x ) 


n 


( 12 ) 

(13) 

(14) 


The form (|T4|) coincides with that assumed in the 3 Pq model. At this 


point one should take into account that the pair in (|I4D is created at any 
point x in the space, and by the derivation this point should lie on the world 
surface of the string, deforming it from the general minimal area shape to 
the minimal area surface passing through the point x, X 4 . 

The probability amplitude for such a deformation is equal to exp(—Ad(i)), 
where AA(x) = crAS is the increase of the area AS. Assuming the latter is 
described by the cone with the height h and radius P, h <C P, one obtains 
an additional factor to be inserted in ( 0 ), namely, 


7T . 


F = exp(— AA) ~ exp(— a—h ) 


(15) 


This factor is similar to the one suggested in the 3 P 0 flux-tube model |2| . 
Integrating over h in (|I4|) yields the final result, which we write in the same 
form as in the 3 Pn model PH. 


H\ SBM) = g 


J d 3 x'ip(x)'ip(x); g = 2T g a/y/n] 7 = 


TgO 


9 _ 

2/i’ 


(16) 


where /i is the constituent quark mass computed through a [|1^]. Taking 
T g ~ (0.2—0.3) fm, a ~ 0.2 GeV 2 and /i ~ 0.35 GeV one obtains 7 ~ 0.3—0.5, 
i.e. the values in the same ballpark as in the phenomenological analysis |3]. Q. 
Several remarks are now in order. Firstly, above the simplified Gaussian 


form of D(x) was used in []T 0 |, which requires the redefinition of T g \ secondly, 
the contribution of magnetic correlators and the term Di in (|lj) has been 
neglected. These additional contributions will be considered elsewhere. 
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4 “Perturbative” pair creation via the hybrid 
or glueball formation 

Two decay mechanisms discussed in previous sections are of purely nonper- 
turbative origin, now we turn to the mechanisms which contain as a limit 
a purely perturbative qq pair creation by a gluon. At large distances one 
has to know how this process is modified by the presence of nonperturbative 
confining fields and to this end we shall use the Background Perturbation 
Theory (BPTh) developed in where the total gluonic held A M is sepa¬ 


rated into valence gluon held a^ and background BA M = B^ + a^. The held 
B M saturates correlator D(x) and contains therefore its own mass scale, while 
perturbation theory is done in powers of ga^. The main physical outcome 
of the analysis of |10| is that the valence gluons are propagating inside the 
him of the string worldsheet, so that all Feynman diagrams in the coordinate 
space can be considered as filled inside by this him on the minimal surface 
with boundaries specihed by quark and gluon trajectories. 

At this point it is clear that one should use the path integral represen¬ 
tation for tha quark and gluon Green’s function, namely the Fock-Feynman- 
Schwinger (FFS) formalism (T!|. The FFS method has proved useful in 
conjunction with the BPTh to study meson, hybrid and glueball Green’s 
functions (see [18] for review). In [2IJ this method has been exploited to cal¬ 


culate mixing between meson, hybrid and glueball states, and in what follows 
we shall pursue the same way to study matrix elements of decays proceeding 
via hybrid and glueball intermediate states. 

We start with the OZI allowed planar pair creation mechanism by a gluon, 
propagating inside the him in a hybrid state ^ H \ anf therefore the matrix 
element for the meson decay via the hybrid states can be written as 

W {H) = Y, X h MH)W n H) (17) 

n 

where is the dimensionless mixing coefficient of the n-th hybrid state 

in the given initial meson, which according to [^Oj can be written as 

\ (MH) = V ™— - (18) 

yj2/i g (n)\Mp - M m \ 

Here M^\ M m are hybrid and meson masses respectively and /i 9 (n) is 
the constituent gluon mass in the n-th hybrids state, computed through the 
string tension as in 


18 , 20 
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The matrix element V^\ introduced in [P0| , is that of the pair creation 
operator Hi, 

Hi = g J g(x, 0 )a(x, 0 )g(x, 0)d 3 x (19) 

between the meson state $^(r) and the hybrid state <f>^j(r 1; r 2 ), where we 
have specihed quark Dirac indices a/3 , 7 S and gluon component index /i (this 
is the 4x4 representation typical for the Bethe-Salpeter wave functions which 
is used to build up 2 j + 1 components of meson and hybrid wave functions, 
for discussion and refs, see H). 


In a simplified form V on can be written as [H 


V on = g J d 3 r</>*^( 0 , r)V(j) ( ' M \r)d 3 r + perm. 


( 20 ) 


where ’’perm” implies another term with gluon emitted by antiquark, and 

r = 7,- 

Similarly the term W^ H) in ([17]) is the decay amplitude of the n-th hybrid 
state into two mesons, which can be written as 


W (H) = Vn, 12 , ^ = r $*(M 1 )( ri )$*(M 2 )( r2)r (H)( ri , r2 ) d 3 fid 3 f2 

feUgWV J 

( 21 ) 

Normalization of wave functions is (summation over repeating Dirac in¬ 
dices, which are omitted, is implied) 


J \$ (H) \ 2 d 3 rid 3 r 2 = 1 , J (<f> (Mi) | 2 d 3 r = 1 (22) 

Finally two-body decay probability is 

Vd 3 h 

w = 2n\WW\ 2 5(Ei -E 2 - M)j— (23) 


As it is seen from S3), (0). the decay probability is strongly enhanced 
when the decaying meson mass is in the vicinity of some hybrid level. One 
should stress that this situation is standard in the mass range above 1.4. 
GeV, where the density of hybrid ground and excited states is fast growing 
with mass, (see e.g. second ref. in ||). 

We now turn to the OZ1 forbidden, i.e. nonplanar qq pair creation via 
valence gluons. The essential step in this mechanism is the creation of the new 








flavour quark-antiquark pair with the new string between them, hence the 
creation operator contains at least a two-gluon exchange. In the confining 
background trajectories of these two gluons are connected by the adjoint 
string (or, equivalently at large N c by a double fundamental string) and 
therefore a new meson is created by the two- or more gluon glueball, which 
is emitted virtually from the original meson. The amplitude for this process 
can be written as follows, 


W (G) 


E 


\T M M2 'W'GMi 
v n v n 


M-E%~ E(M 2 ) 


where notation is used 


V MM, = )| A (G)|0(«)) , 

V GM, = 


(24) 


(25) 


and the two-gluon glueball vertex is 

A i% 0 i,x 2 ) = K { Su2( x ^ x 2)(lu 1 S {f \x l ,x 2 )^ U2 ) a p 

2 ) = (l^S {9 \y 1 ,y 2 )'y l , 2 ) a (3^ ltl2 (yi,y2) (26) 


Here TV) is the glueball wave function, and SV/) _ the Green’s function of 
quarks with flavor (g, /). 

When all distances \yi — Xk\, \x\ — x 2 \, Ui — y-i are small as compared to 
T g , the two-gluon Green’s function reduces to the two-gluon exchange of free 
gluons. This limiting perturbative mechanism is known for a long time |2T| . 
In the opposite limit only the lowest mass term survives in the sum ([h|). 

The most interesting case occurs when the denominator in (^4|) becomes 
small, which is possible when the mass of some glueball state is close to the 
mass of the emitted meson M\. This amplification may thus occur in the 
0 ++ channel for M\ around 1.5 GeV, or in other channels for Mi > 2 GeV 
22j . A similar mechanism may take place in hadron-hadron scattering with 
creation of the cc states in the mass range 3-4 GeV, where also glueball 


states are predicted in lattice and in analytic calculations [22 


5 Conclusion 

In the present paper the first qualitative step was done aimed at the system¬ 
atic derivation of strong decay amplitudes from the QCD Lagrangian. We 
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have considered three decay mechanisms, where nonpertnrbative contribu¬ 
tion is very important. 

The first one, the CDM, was already successfully checked in heavy-light 
meson decays . The next step would be to apply CDM to the pionic and 


kaonic decays of light-light mesons, and investigate double and more-pion 
decays, which are given by the nonlinear Lagrangian (|8|). The second mech¬ 
anism, the SBM, turns out to be mostly 3 P 0 (additional terms due to the 
smaller nonscalar component in the quark Green’s function S considered in 
|TT| have been neglected above) and its amplitude is close to the phenomeno¬ 
logical fits [|]. The third type of mechanisms with hybrid or glueball in the 
intermediate state is the nonperturbative background generalization of the 
original purely perturbative mechanisms of the 3 S i type 
gluon exchange j2l] , 


and of two- 


respectively. It is argued that a strong enhancement of 
decay amplitudes is possible when the corresponding levels of a hybrid or a 
glueball are close to the mass of the original or the final meson respectively. 
The paper does not contain quantitative predictions, which are planned for 
the future, but rather is concentrated on the general discussion of possible 
decay mechanisms as they emerge from the basic QCD Lagrangian. 

One of immediate extensions of the present result is the inclusion of the 
baryon decays, where all three mechanisms discussed above are present in 
the same form with the only replacement of the simple mesonic string by the 
Y -shaped baryonic string. 

Finally, one should have in mind that all three mechanisms discussed 
above enter the decay amplitude additively and therefore one can expect 
in general some interference effects, which make the analysis of data more 
complicated and perhaps more interesting. 

The author is grateful to Yu.S. Kalashnikova and V.I. Shevchenko for 
valuable discussions. The financial support of INTAS grants 00-110 and 
00-366 and RFBR grants 00-02-17836 and 00-15-96786 is gratefully acknowl¬ 
edged. 
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